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Abstract 



■ New heavy charged lepton production and decay signatures at future 

' electron-positron colliders are investigated at ,/s — 500 GeV. The conse- 

quences of model dependence for vector singlets and vector doublets are 
C**) ■ studied. Distributions are calculated including hadronization effects and 

' experimental cuts that suppress the standard model background. The 

pH , final state leptonic energy distributions are shown to give a very clear 

signature for heavy charged leptons. 

2^ ■ 1 Introduction 



The possibility of new heavy charged leptons is present in many extentions of 
the Standard Model (SM). In grand unified theories beyond SU(5) we have 
new fermionic fields pQ. In superstring theories with an Eq gauge group after 
compactification, the fermions in the fundamental 27 representation have new 
leptonic doublets and new neutral lepton singlets [3> . In left-right symmetric 
models 0] we have also new mirror fermions that are expected to be much 
heavier than the presently known fermions from the SM. Extended electroweak 
gauge groups with three anomaly free families 5 contain new fermions. We can 
also have a new fourth family, similar to the electron-muon-tau families, but 
with much higher mass states. 

It is well known that experimentally no clear signal of such fermions was 
found so far If these particles exist, the latest experimental limits imply 
that they must have a mass greater than 100 GeV. However, the possibility 
of generalizing the SM into a deeper model is a strong motivation for both 
experimental searches at the next colliders as well as for theoretical studies. 
The recent experimental confirmation of light neutrino masses and oscillations 
\7\ also raises the question of new possible heavy neutral leptons. These neutral 
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heavy leptons could be responsible for the see-saw mechanism for small neutrino 
masses. As the presently known leptons appear in neutral and charged states, 
it is reasonable to suppose that new heavy states will show the same behavior. 
The new heavy charged leptons could be searched for in future high energy 
hadron-hadron colliders 0, lepton-hadron colliders 9,, lepton-lepton colliders 
[TU] and photon-electron JT] colliders. In this paper we will concentrate our 
attention on the electron-positron option at y/s = 500 GeV. This will be the 
first energy option for the linear collider at NLC ^2 as well as for the TESLA 
[T3") project at DESY before they move to higher center of mass energies and 
luminosities. Most of our results can be scaled for higher energy options. We 
have considered an yearly integrated luminosity of C = 10 5 pb _1 . 

The new heavy charged and neutral leptons have been studied intensively 
by many authors PElEHHElIHICllHIini EDI CU C3I CS1 IIU For heavy lepton 
masses up to Ml = y/s/ 2 the total cross section is dominated by pair production 
due to photon exchange. So this mechanism is quite model independent. We 
can extend the heavy lepton mass range to y/s > Ml > y/s/ 2 by looking for 
single heavy lepton production. In this case, and also in the decay of heavy 
leptons, there is a strong model dependence on the mixing between light and 
heavy leptons as well as in the heavy-to-heavy lepton couplings. 

In the present work we study the dynamical and kinematical properties for 
production and also for the decay of new heavy charged leptons into presently 
known particles. These distributions are compared with the expected back- 
ground coming from the SM and we show how the signal can be clearly separated 
from the background with a great significance, independent of hadronization ef- 
fects. Another point discussed is the model dependence for new heavy charged 
lepton experimental search. 

This paper is organized as follows. In section 2 we discuss the models adopted 
and their parameters. In section 3 we present the results from the various pro- 
cesses involving heavy leptons, studied in this paper, comparing the distribu- 
tions from SM, vector singlet model (VSM) and vector doublet model (VDM). 
Finally, conclusions are shown in section 4. 

2 The models 

The high precision data on the SM parameters imply strong limits on new heavy 
gauge bosons. So, our first hypothesis is that in the next e^e + collider at y/s = 
500 GeV the mechanism for production and decay of new heavy charged leptons 
will be dominated by the presently known electroweak gauge bosons. The second 
hypothesis is that the new heavy leptons can mix with the light leptons in each 
family. The mixing angle parameters are constrained by the Z properties and by 
asymmetries in lepton-lepton scattering. They imply a small mixing j!4j . of the 
order of sin 2 9 m ix < 0.005 with 95% c.l. Throughout this paper we will consider 
all mixing angles (#f, fl and 9 V L R ) between the new and the ordinary leptons 
(charged and neutral) equal to the present upper bound sin 2 9 mix = 0.005. The 
next theoretical point concerns the quantum numbers that can be assigned to 
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Table 1: Vector singlet model couplings between bosons and fermions with N 
as a Dirac field 
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new heavy leptons. Here we have a number of possibilities, but we will restrict 
our analysis to only two models: new vector singlets and new vector doublets. 
For the VSM the new charged heavy lepton (written as L) is a singlet, as well 
as the new neutral lepton (written as N). For the VDM both new heavy leptons 
are placed in doublets. There are other possible models for new heavy leptons 
such as mirror and fourth generation models. We have restricted our analysis 
to the VDM and VSM since they give the more important features of model 
dependence with fewer parameters than other models. 

The interactions are then given by the Lagrangians for neutral and charged 
currents: 



ten = ~^ Z ^ {4 ~ 5a7 5 ) $j + h.c. (1) 

Ccc = -^ W ^ (dv - dhn) $i + h.c. (2) 

where i and j are the possible combinations of light and heavy leptons. The 
couplings depend on mixing parameters that are given in table 1 for the VSM, 
and in table 2 for the VDM. In the next section we will discuss the role of heavy 
neutrinos in the heavy charged lepton decays. In tables 1 and 2, we present 
the couplings for the case of new heavy neutrino states of the Dirac type. In 
this case a global family symmetry is present and no lepton number violation 
occurs. 



The heavy neutrinos could also be of the Majorana type. Then there are 
small changes in their couplings due to the identity ^n^^^n = 0, i.e., pure 
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Table 2: Vector doublet model couplings between bosons and fermions with N 
as a Dirac field 
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vector couplings are zero. We have in each family the same couplings of the 
heavy Majorana neutrino with ordinary lcptons. In the Majorana case, we 
have the possibility that heavy charged lepton decays could also violate lepton 
number conservation |15| due to the chain, 



NW 
I 

— >l± W T 



(3) 



where £,i = e,n,T, so that the cross sections are the same for dilepton final 
states of the type e~e+ — > e ± + + W ± + W T . 

In all calculations presented in this paper we have considered TV as a heavy 
Majorana neutrino that couples to all families, with mass Mn — 100 GeV. 
We are assuming the heavy neutral lepton mass to be smaller than the heavy 
charged lepton mass as in the light lepton mass spectrum for the three known 
lepton families. 

For simplicity we have concentrated on the first (electron) family and our 
results can be generalized to the other families. 



3 Results 

In this section we present the main channels, cross sections and distributions, 
comparing the signal from VSM and VDM and the SM background. This pro- 
cedure allow us to obtain a clear signature for heavy charged leptons. As many 
of these processes involve a large number of diagrams, we have employed the 
CalcHep ^7] package. The VSM and VDM were implemented in the unitary 
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and Feynman gauges and both results were carefully checked and compared. 
So, the results here presented are complete first order calculations. 

In Figs, la and lb we give the total and partial widths for heavy charged 
lepton decay as a function of its mass. The enhancement for the VDM model, 
shown in Fig. la, comes from the heavy-to- heavy transition in the channel 
L e — > NW as shown in Fig lb and it starts when = M/v + Mw- This 
channel is suppressed in the VSM due to the sin 2 6 m i X factor that is present in 
the W — L — N vertex, as shown in the last line of table 1. For the VDM, the 
corresponding term is maximum for this vertex, as shown in table 2, strongly 
suppressing the decays L e — > e + Z and L e — > v e + W. Of course, if the neutral 
heavy lepton mass increases, this effect is smaller. We can have important model 
dependence for the signatures of charged heavy lepton decays. 

For pair production of heavy charged leptons we have three diagrams, but 
the result is dominated by photon exchange. This is the leading process for 
masses up to \fsj2. The single heavy lepton production gives a significative cross 
section for masses > \fsf1 an d depends strongly on the mixing parameters. The 
corresponding total cross sections for real pair and single charged heavy lepton 
production are given in Fig. 2. Both the VDM and the VSM give the same 
results. In order to improve our calculations as close as possible to experimental 
situations, we have employed the following detector cuts for final state particles: 
a minimum 5 GeV energy cut for all final leptons; an angular cut \cosO\ > 0.995 
for all particles, except for the final gauge bosons, relative to the initial electron 
beam direction and an invariant electron-positron mass cut > 5 GeV. 

For single heavy lepton production one can easily compute cross sections 
for the three particle final states eeZ and veW by multiplying the results of 
Fig. 2 by the corresponding branching fractions. As we are mainly interested 
in comparing different signatures with the SM background, we have performed 
the full first order calculation. The Feynman diagrams for the signal are given 
in Figs. 3a and 3b. We notice that in Fig. 3a, each diagram appears twice. For 
the SM background we have sixteen diagrams for the first channel and twelve 
diagrams for the second. This procedure of including all diagrams, allows us to 
have a reliable result not only for cross sections but also for final state particle 
distributions [TSJ. These large number of diagrams can be taken into account by 
employing high energy algebraic programs, such as CalcHep [T?j- The resulting 
cross sections as a function of Ml for the processes + e + — ► + e T + Z 
and e + + — ► e + (e~) + u e (i> e ) + W~(W + ) are shown in Figs 4a and 4b. 
The curves labeled VDM and VSM correspond to the signal terms only. The 
first channel has only the charged heavy lepton signal (with both charges) but 
the second channel has both neutral and charged heavy leptons (with only one 
charge) contributions, as shown by Figs. 3a and 3b. We show in Fig. 4b 
the role of the heavy neutrino contribution to the signal total cross section, by 
computing separately the heavy charged lepton contribution only. Since the 
dominant process is the real production of heavy leptons, our results are gauge 
invariant. These are the main signatures for a single heavy charged lepton with 
mass in the 250 — 500 GeV range. 

For pair production with masses in the 100 — 250 GeV range we must look 
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Figure 1: Total heavy lepton width (la) and partial widths (lb) for VSM and 
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Figure 2: Total cross section for single and pair production of heavy leptons in 
e- + e+ at V« = 500 GcV (sin 2 9 mix = 0.0052) for VSM and VDM. 
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Figure 3a 




Figure 3: General Feynman graphs for charged heavy lepton contribution in 
e~ + e+ — ► e ± + e T + Z (3a) (8 dominant diagrams) and e~ + e + — ► e+ + 
v e + W~ (3b) (6 dominant diagrams) with Mn = 100 GeV. 
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Figure 4: Total cross sections for e + e + — ► e ± + e T + Z (4a) and e + e + — ► 
e+(e") + i/ e (P e ) + W-(W+) (4b) at ^/i = 500 GeV with M N = 100 GeV. 
In fig. 4b we have separated the pure9signal N + L e from the pure signal L e 
contributions. 
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Figure 5: Total cross sections for e + e + — > e ± + e T + Z + Z and e + + e — > 
e+(e") + v e {v e ) + W~(W+) + Z at y/s = 500 GeV with M w = 100 GeV. 

for final states with four particles. For example in Fig. 5 we show the cross 
sections as a function of Ml for the channels e+(e~) + v e {v e ) + W~(W + ) + Z 
ande ± + e =F + Z + Z. The curves labeled VSM and VDM, include only the signal 
from heavy neutral and charged leptons. We notice that the signal is above the 
SM. For the VDM model and Mm— 100 GeV, the contribution of new heavy 
leptons to the total cross section is the same as in the VSM up to the region 
ML e = Mm + Mw- Then the dominant decay mode is L e — > N + W and the 
other channels become suppressed. 

In the experimental search for heavy charged leptons a fundamental point is 
the choice of distributions that can enhance the signal from the SM background. 
A simple distribution is the final state invariant mass M e z, as shown in Fig. 6 
for the e~e + Z channel, for the complete VSM extension (signal + background). 
This distribution requires the experimental reconstruction of the Z hadronic 
jets. There are intrinsic difficulties and uncertainties in the experimental and 
theoretical reconstruction of hadronic events. 

In order to account for the hadronic effects in our study we have employed 
the package PYTHIA 6.1 ^B] to generate the quark hadronization and decay 
from real W and Z production. 

An important point that we call attention in this paper is the kinematical 
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Figure 6: Reconstruction of the invariant mass M e z for e+ + v e + W + Z at 
Vs = 500 GeV with M N = 100 GeV for the complete VSM extension (signal + 
background) . 
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Figure 7: Final lepton energy and invariant mass distributions in the e + + v e + 
W~ channel at = 500 GeV, for M N = 100 GeV and M La = 450 GeV. 



relation for a primary electron (positron) produced with a L e that must satisfy: 

Ee ±AE=^^n (4) 

For a sharp Ti e we will have a sharp energy distribution of the primary 
lepton produced together with the new heavy charged lepton. This property 
suggests interesting signatures for heavy charged lepton searches by looking 
at final state charged lepton distributions. This is far more effective than the 
reconstruction of invariant masses or hadron jets from W and Z hadronic decays. 
Analyzing the Fig. 7, for v e + e + + W~ channel, we can see that the peak 47.5 
GeV (240 GeV) at the final electron (neutrino) energy distribution is related to 
the 450 GeV (100 GeV) heavy charged (neutral) mass. 

We also have reconstructed the final state electron energy distribution in the 
e ± +e =F + hadrons channel. The background and Z hadronization were generated 
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with the PYTHIA package as mentioned above. A clear separation between 
signal and SM background is given by 2-dimensional plots instead of usual 1- 
dimensional plots. We have considered the 2-dimensional phase space plot of 
the electron and positron energy in the e +e^+ hadrons channel. According to 
equation 4, the electron and positron energy of the events coming from the signal 
are peaked at 160 GeV for a 300 GeV charged lepton mass. This correspond 
to the two crossed lines in Fig. 8a since the process is symmetric. The other 
concentration at the inclined line is from the standard model Z — > e + + e~ 
decay. For Ml c = 450 GeV the result is given in Fig. 8b with the dense lines 
of the signal shifted to the corners of the kinematical region. For the e + + v e + 
hadrons channel the analogous plot is given by the final neutrino missing energy 
versus the final positron energy. For this channel we have the contribution of 
both the neutral and charged heavy leptons, as shown in the Feynman diagrams 
of Fig. 3b. In this case we have one line in the missing energy that comes from a 
neutral heavy lepton decay {Mm = 100 GeV) and the other line in the electron 
energy comes from the charged heavy lepton decay (Ml c = 450 GeV). The 
result is given in Fig. 9 which is a detail of the allowed phase space. We call 
attention to the fact that the heavy neutral lepton line is more defined than the 
heavy charged lepton line. This is due to the higher cross section for the neutral 
heavy lepton, as shown in Fig. 4b and to the different widths. For M Lc = 300 
GeV we have a similar figure with an horizontal line at the same position since 
we maintain Mm = 100 GeV and a vertical line at E e — 160 GeV. 

Besides the general cuts mentioned before, the phase space figures suggest 
that one can perform additional more restrictive cuts in order to enhance the 
signal over the background. The numbers of events were obtained for an integral 
luminosity of L = 10 5 pb _1 and an upper bound for the mixing angle sin 2 9 m ix — 
0.005. Table 3 gives the total number of events for the signal s (VSM) and for 
the SM background 6, as well as the statistical significance s/Vb after all cuts. 
As the planned luminosities could be increased ^2] by a factor of 3, one could 
explore even lowers values for the mixing angle. 



For the e + e =F + hadrons channel, the resulting number of events after all 
cuts is given in Fig. 10a for an integrated luminosity of 10 5 pb _1 . The analogous 
result for the e + + v e + hadrons channel is given in Fig. 10b. 

The lepton number violating decay of heavy charged leptons is significant 
only in VDM. If Ml c > Mm + Mw and Mm > Mw we can have processes of 
the type e + + e~ — > e~ + hadrons due to the chain: 



NW + 
I 

— /*- W+ (5) 
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Figure 8: Positron energy versus electron energy event plot in the e ± + 
hadrons channel for M^ e = 300 GeV (8a) , and the same distribution for Ml £ = 
450 GeV (8b) &tJa = 500 GeV for VSM. 



e - + e * ^ e * + v + W - 




Figure 9: Signal region of interest of the phase space for electron energy versus 
missing energy event plot in the e + + v e + hadrons channel for Ml c — 450 GeV 
at V« = 500 GeV with M N = 100 GeV for VSM. 
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Figure 10: Number of events versus final charged lepton energy in the e ± + 
e T + hadrons (10a) and in e+ + v e + hadrons (10b) at y/s = 500 GeV with 
M N = 100 GeV for VSM. 
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Table 3: Number of events after hadronization and selection cuts for signal 
(s)(VSM) + background (b) (SM), background, significance and energy cuts. 
The cut E e + + E e - (i?) > 251 GeV eliminates the inclined line in phase space 
corresponding to Z — > e + + e~ (W — > e + v). 



e + + e — ► e + + e + hadrons 


Lepton masses 


s + & 


b 


a/Vb 


Cuts 


Reconstructed mass 


M Le = 300 GeV 


764 


251 


30 


E e+ + E e - > 251 GeV 
158 < E e ± < 162 GeV 




M Le = 450 GeV 


396 


205 


13 


E e+ + E e - > 251 GeV 
45 < E e ± < 49 GeV 


M Lc 


e + + e — ► e + + u e + hadrons 


M Lc = 300 GeV 
M N = 100 GeV 


1357 


599 


31 


E e+ + > 251 GeV 
158 < E e + < 162 GeV 


M Ls 


M Lf = 450 GeV 
M N = 100 GeV 


477 


218 


18 


E e + + $ > 251 GeV 
45 < E e+ < 49 GeV 


M Lc 


M Le = 300 GeV 
M N = 100 GeV 


2582 


599 


81 


E e + + $ > 251 GeV 
239 < $ < 241 GeV 


M N 


M Le = 450 GeV 
M N = 100 GeV 


2515 


599 


78 


E e + + > 251 GeV 
239 < # < 241 GeV 


M N 



The total cross section is given in Fig. 11. Here again the single final 
state electron energy distribution can give information on the mass of the heavy 
charged lepton as shown in Fig. 12, but with a much larger r^e when compared 
with the VSM. The reconstructed invariant mass M^ww gives the same infor- 
mation, but requires the hadronic reconstruction of the invariant masses and/or 
jets of two Ws with higher uncertainties. 

4 Conclusions 

In this paper we have investigated different signatures for heavy charged lep- 
ton detection at the next electron-positron colliders. We have presented cross 
sections and distributions for the ^fs = 500 GeV option for NLC and TESLA 
colliders, but our conclusions are very similar for higher energies. For vector 
charged leptons we notice that vector doublets and singlets give the same results 
if a new heavy neutrino in the same representation has Mjy — Ml c ■ However, 
if the heavy neutrino mass is smaller than the new heavy charged lepton, then 
vector doublets and singlets will show a very different behavior. The decay 
L e — > N + W could be dominant for vector doublets and suppressed for vec- 
tor singlets. If the new heavy neutrino is of the Majorana type, we will have 
the possibility of lepton number violation in the heavy charged lepton decay. 
This property allows the experimental determination of the Dirac or Majorana 
nature of heavy neutral leptons. 

For new heavy charged leptons with mass up to M^ e w y/s/2, the dominant 
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Figure 12: Electron energy distribution in the e + \l + W + + W + channel at 
= 500 GeV with M N = 100 GeV for VDM. 
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process is pair production via photon exchange. Then the final lcpton energy 
distribution in the processes e ± + e T + Z and e+(e~) + v e (v e ) + W~(W + ) 
allows a clear separation from the SM background. If the new heavy charged 
lepton has a mass in the range y/s/2 < Ml c < \A then single production is a 
feasible process. In this case we call attention to the kinematical properties of 
the associated primary lepton produced with the new charged lepton, according 
to the equation 4. Its energy distribution can allow a very precise determination 
of the new heavy lepton mass without the looses and errors that came from the 
reconstruction of hadronic events from W and Z decays. This kind of analysis, 
in the future e~ + e~ colliders, can lead to much stronger bounds on new heavy 
lepton parameters if no signal is found . 
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